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Abstract 

The measured particle ratios in central heavy-ion collisions at RHIC-BNL are investigated within 
a chemical and thermal equilibrium chiral SU (3) ct — lo approach. The commonly adopted noninter- 
acting gas calculations yield temperatures close to or above the critical temperature for the chiral 
phase transition, but without taking into account any interactions. Contrary, the chiral SU{'i) 
model predicts temperature and density dependent effective hadron masses and effective chemical 
potentials in the medium and a transition to a chirally restored phase at high temperatures or 
chemical potentials. Three different parametrizations of the model, which show different types of 
phase transition behaviour, are investigated. We show that if a chiral phase transition occured in 
those collisions, "freezing" of the relative hadron abundances in the symmetric phase is excluded 
by the data. Therefore, either very rapid chemical equilibration must occur in the broken phase, or 
the measured hadron ratios are the outcome of the dynamical symmetry breaking. Furthermore, 
the extracted chemical freeze-out parameters differ considerably from those obtained in simple 
noninteracting gas calculations. In particular, the three models yield up to 35 MeV lower temper- 
atures than the free gas approximation. The in-medium masses turn out differ up to 150 MeV 
from their vacuum values. 



INTRODUCTION 



Thermodynamical equilibrium calculations of particle production in high energy particle- 
and nuclear collisions have been carried out for a long time |jl|, |, |, |], |^ ||, ^, |10|, |lT 



Recently hadron abundances and particle ratios have been measured in heavy-ion collisions 
from SIS, AGS, SPS to RHIC energies. These data have revived the interest in the extrac- 
tion of temperatures and chemical potentials from thermal equilibrium "chemical" model 
analyses. The experimentally determined hadron ratios can be fitted well with straightfor- 
ward noninteracting gas model calculations 0, |^, ^ |TD|, [TI|, |T3] , if a sudden breakup of 
a thermalized source is assumed and once the subsequent feeding of the various channels by 
the strongly decaying resonances is taken into account. From the freeze-out fits one has 



constructed a quite narrow band of freeze-out values in the T — fiB plane (see e.g.lTJ, p^). 
The extracted freeze-out parameters are fairly close to the phase transition curve for SPS 
and RHIC energies. However, when we are indeed so close to the phase transition or to 
a crossover as suggested by the data for T and yU^, we can not afford to neglect the very 
in-medium effects we are after - and which, after all, do produce the phase transition. Thus, 
since noninteracting gas models neglect any kind of possible in-medium modifications they 
can not yield information about the phase transition. 

Therefore, we will employ below a relativistic selfconsistent chiral model of hadrons and 
hadron matter developed in |I6|. This model can be used as a thermodynamically 

consistent effective theory or as a toy model, which embodies the restoration of chiral sym- 
metry at high temperatures or densities. Therefore the model predicts temperature and 
density dependent hadronic masses and effective chemical potentials, which have already 
been proposed and considered in P, [T^, Thus, using the chiral SU{3) model 

we can investigate, whether the freeze-out in fact takes place close to the phase transition 
boundary (if it exists) and if the extracted T,fiB parameters are strongly model dependent. 
Depending on the chosen parameters and degrees of freedom different scenarios for the chiral 
phase change are predicted by the model: Strong or weak first order phase transition or a 



crossover. The transitions take place around = 155 MeV |T^, |21|, which is in qualita- 
tive agreement with lattice predictions [E^ for the critical temperature for the onset of a 



deconfined phase which coincides with that of a chirally restored phase [23 
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MODEL DESCRIPTION 



The chiral SU{3) model is presented in detail in |T^, |T6|. We will briefly introduce the 
model here: We consider a relativistic field theoretical model of baryons and mesons built 
on chiral symmetry and broken scale invariance. The general form of the Lagrangean looks 
as follows: 

^ = -^kin + -^BW + -^VP + -^vcc + ^0 + ^SB- (l) 

W=X,Yy,A,u 

£km is the kinetic energy term, £bw includes the interaction terms of the different baryons 



with the various spin-0 and spin-1 mesons (see |T6[ for details). The baryon masses are gen- 
erated by both, the nonstrange a {< qq >) and the strange ( {< ss >) scalar condensate. 
£vp contains the interaction terms of vector mesons with pseudoscalar mesons. £vec gener- 
ates the masses of the spin-1 mesons through interactions with spin-0 fields, and Co gives the 
meson-meson interaction terms which induce the spontaneous breaking of chiral symmetry. 
It also includes a scale-invariance breaking logarithmic potential. Finally, C^b introduces an 
explicit symmetry breaking of the U{1)a, the SU(3)y, and the chiral symmetry. All these 
terms have been discussed in detail in |T6|. 

The hadronic matter properties at finite density and temperature are studied in the mean- 
field approximation [^. Then the Lagrangean (P becomes 



i 

Cvec = -m^^u +-m.^(l) +g^[uj + 2(f) ) 



Vo = lkox\cr' + e) - h{a' + C^)^ - + - hxa\ (3) 



\x^ In ^ - ^-x^ In ^ 
4 Xo 3 (T^Co 



+ k,x' + 7r In ^ - -x' In ^ (4) 



VsB = (—) mlf^a + {\f2m\fK - -^^If^ 



)C 



where is the effective mass of the hadron species i. a and ( correspond to the scalar 
condensates, u and represent the non-strange and the strange vector field respectively, 
and X is the scalar-isoscalar dilaton field, which mimics the effects of the gluon condensate 
25| . Only the scalar (Cbx) and the vector meson terms (Cbv) contribute to the baryon- 



meson interaction, since for all other mesons the expectation value vanishes in the mean-field 
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approximation. The grand canonical potential Q per volume V as a function of chemical 
potential fi and temperature T can be written as: 

^ = -£.ec - A - CsB - V.ac T T ^ ^ I d'k [lu (l ± e" ^ [^*(^)-^*] ) ] , (5) 

with the baryons (top sign) and mesons (bottom sign). The vacuum energy V^ac (the po- 
tential at pb = 0,T = 0) has been subtracted in order to get a vanishing vacuum energy. 
7i denote the hadronic spin-isospin degeneracy factors. The single particle energies are 
^ii^) = V^i + "^i^ the effective chemical potentials read fi* = pi — Qi^jUJ — gcj^icj). 
The mesonic fields are determined by extremizing ^(/x, T = 0). The density of particle i can 
be calculated by differentiating Vt with respect to the corresponding chemical potential /i,. 
This yields: 



Pi = li 



(6) 



All other thermodynamic quantities can also be obtained from the grand canonical potential. 
In the present calculation the lowest lying baryonic octet and decuplet and the lowest lying 
mesonic nonets are coupled to the relativistic mean fields. Depending on the coupling of the 
baryon resonances (the decuplet) to the field equations, the model shows a first order phase 
transition or a crossover (for details see [^). We will use three different parameter sets: 
Parameter set CI treats the members of the baryon decuplet as free particles, which yields 
a crossover behaviour. Parameter sets CII and GUI include also the (anti)-baryon decuplet 
as sources for the meson field equations. They differ by an additional explicit symmetry 
breaking for the baryon resonances along the hypercharge direction, as described in 



for the baryon octet. This is included in CII and not used in CIII. This leads to a weak 
first order phase transition at yU = for CII and two first order phase transitions for CIII, 
which can be viewed as one strong first order phase transition. Heavier resonances up to 
m = 2 GeV are always included as free particles. The resulting baryon masses for CI and 
CIII are shown in fig. |I|. We observe a continous decrease of the baryon masses for CI starting 
at T ^ 150 MeV. In contrast, CIII shows two jumps around T = 155 MeV. The critical 
energy densities, the entropy densities and the transition temperatures for /i^ = /i^ = 
{fj,q = fj,B/'i,fJ's = — Ps) are specified in table |. 
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FIG. 1: Baryon octet masses as function of temperature for vanishing chemical potential. Left CI, 
right cm. Note the continuous change of the masses starting around T = 150 MeV. In contrast 
for cm two phase-transitions occur around ~ 155 MeV. These result from the separate jumps 
in the nonstrange (a) and the strange (C) condensate. 





e"/eo 


e+/eo 


s [fm ^] 


s+[/m ^] 


TjMeV] 


CII 


2.8 


7.2 


2.8 


6.7 


156.3 


cm - 1st PT 


2.3 


8.3 


2.4 


7.9 


153.4 


cm - 2nd PT 


10.5 


17.1 


9.8 


15.7 


155.5 



TABLE I: Energy density, entropy density and phase transition temperatures for CII,CIII, fig = 
fj,s = 0. The (— ), (+) signs refer to an approach to the phase transition from below and above, 
respectively. Tc denotes the phase transition temperature, eo = 138.45 MeV/fm^ denotes the 
energy density of nuclear matter in the ground state. 

PARTICLE RATIOS IN THE CHIRAL SU{3) x 5C/(3) MODEL 

Since the chiral SU(3) model predicts density and temperature dependent hadronic masses 
and effective potentials, in contrast to noninteracting models, the resulting particle ratios 
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and therefore the deduced freeze-out temperatures and baryon chemical potentials are ex- 
pected to change Hence in the following, we identify combinations of temperatures and 
chemical potentials that fit the observed particle ratios in the chiral model. In all calcula- 
tions the value of the strange chemical potential /^^ is chosen such that the net strangeness 
fs = 0. We are looking for minima of with 

i * 

Here r*^^^ is the experimental ratio, r^"'^'^'- is the ratio calculated in the model and cTj repre- 
sents the error in the experimental data points. We use the same ratios as in |Q: p/p, A/A, 
S/H, n-/n+, K-/K+, K-frr-, K*/h-, K*/h-. 

Even though the only parameters in a thermal and chemical equilibrium approach on first 
sight are the temperature and the baryon chemical potential, there exist further unknowns: 
On the one hand, some decays of high mass resonances are not well known and on the 
other hand the effect of weak decays in the experiments strongly depends on the detector 
geometry and on the reconstruction efficiency of the experiments. The feeding correction 
from the strong and electromagnetic decays of the hadronic resonances used here employs 



the procedure used in the UrQMD model Weak decays are not considered here. 

We rather focus on the principal question whether an interacting chiral SU(3) approach 
with m* 7^ rrivac can at all describe the particle yields at RHIC. Fine tuning of the by 
adjustment of the weak decay scheme is not our intention. Even though it has been shown 
p9| that values may be improved by including weak decays. 

To compare the quality of the fits obtained in the chiral model with those obtained from 
the noninteracting gas approach, we set all masses and chemical potentials contained in the 
chiral model to their vacuum values and again use the same UrQMD feeding procedure as for 
the interacting model. This yields the ideal gas denoted igpFM- We find that the resulting 
ideal gas ratios are not identical but comparable to those obtained in the literature |P, |2^, 



29| , [30|] . The differences should only result from a different treatment of weak interactions 



and from the uncertainty in the decay scheme of high mass resonances. 
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RESULTS FOR AU + AU COLLISIONS AT RHIC 



First, we find that a reasonable fit of the measured particle ratios at RHIC is possible in 
all three phase transition scenarios of the chiral model and the ideal gas case with comparable 
quality. 

Second, the resulting freeze-out values depend on the model employed, i.e. crossover, 
weak first order, strong first order or free thermal gas. 

Third, a reasonable description of the data is impossible above Tc in the models showing 
a first order phase transition. This shows that no direct freeze-out from the restored phase 
is observed. 

Figure |] shows the value of in the T — fiB plane for the crossover case and for the 
strong first order phase transition. We see that the best fit T — /i^ values differ in both 
models. Furthermore, in the crossover case is well behaved as a function of T and /i^. In 
contrast, the model with a strong first order phase transition shows a very steep increase of 

at the phase transition boundary: the quality of the fit decreases drastically due to the 
jump of the effective masses at the phase transition boundary. Above the values are 
inacceptable, > 500. 




MB[MeV] /iB[MeV] 



FIG. 2: contours in the T — plane for CI (left) and CIII (right). Data are taken from 
On the left, the adiabatic path (constant entropy per net baryon S/A), corresponding to expanisve 
cooling of an ideal fluid, is also shown, /ug is chosen such that = 0. 
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937 


931 
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164 


177 


142 




P/P 


0.65(7) [STAR], 0.64(8) [PHENIX] 
0.60(7) [PHOBOS], 0.61(6) [BRAHMS] 


0.640 


0.648 


0.652 


0.629 


0.629 


A/A 


0.77(7) [STAR] 


0.714 


0.695 


0.702 


0.721 


0.753 




0.82(8) [STAR] 


0.787 


0.731 


0.743 


0.834 


0.894 




1.00(2) [PHOBOS], 0.95(6) [BRAHMS] 


1.000 


1.000 


1.000 


1.000 


1.007 


K-/K+ 


0.88(5) [STAR], 0.78(13) [PHENIX] 
0.91(9) [PHOBOS], 0.89(7) [BRAHMS] 


0.919 


0.914 


0.915 


0.916 


0.894 


K-/TT- 


0.15(2) [STAR] 


0.183 


0.168 


0.168 


0.179 


0.145 


p/vr" 


0.08(1) [STAR] 


0.082 


0.084 


0.078 


0.083 


0.078 


K*/h- 


0.058(17) [STAR] 


0.055 


0.049 


0.049 


0.046 


0.032 


K*o/h- 


0.060(17) [STAR] 


0.049 


0.044 


0.044 


0.041 


0.037 



TABLE II: Chiral fit of the particle ratios measured at RHIC at y/s = 130 GeV. Data and BMRS-fit 
taken from ]|6|. 

The resulting best-fit particle ratios, x^-values and thermodynamic quantities are shown 
in table ^ and figure ^ 

The values for the chiral model are: Xci ~ 5.50, Xcii ~ 5-73 and Xciii ~ 5.40. Thus, 
all three parameter sets describe the data equally well. Furthermore, the agreement is as 
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good as in the noninteracting gas calculation {x% = 5.72 0, Xffm = 5.66). The best fit 
T — fiB parameters vary quite considerably between the different models. The noninteracting 
gas calculation yields T = 187.6 MeV and /i^ = 44.1 MeV. These freeze-out values can be 
compared to those obtained in other ideal gas calculations: T = 174 MeV, /i^ = 46 MeV 
in §, T = 165 MeV, /is = 41 MeV in |Sg and T = 190 MeV,/iB = 45 MeV in 0. The 
crossover case in the interacting chiral model (CI) yields T = 170.8 MeV, /is = 48.3 MeV. 
Very strong deviations are found for the models with a first order phase transition (CII,CIII): 
The freeze-out temperatures are T = 155 MeV (CII) and T = 153.3 MeV (CIII), more than 
30 MeV lower than for igpFM- The fitted baryon chemical potentials fiB increase by about 
7—10 MeV. These T — fiB pairs are very close to the phase boundary (CII) or even right on 
it (CIII) and are about 10 MeV higher than the values obtained at SPS-energies ||3^. Mainly 
due to the different freeze-out temperatures the values of the corresponding thermodynamic 
quantities vary between the different approaches. However, the energy per particle E/A 
is approximately 1 GeV in all cases. This 'unified freeze-out condition' has already been 



proposed in 
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FIG. 3: Particle ratios calculated with CI (left) and CIII (right) compared to RHIC data as 
compiled in [^. 



The fact that the freeze-out appears right at the phase boundary or at crossover implies 
that there are large in-medium corrections, in particular for the effective masses, a phe- 
nomenon observed already in ||T^ . The effective masses shown in figure ^ are shifted up to 



15% from their vacuum values. However, all the interacting models show similar values for 
the effective mass of a given hadron. The strongest in-medium modifications are observed 
for the nonstrange baryons (Am*/mj ^ 15%). Mesons and strange baryons show smaller 
changes of the effective masses, e.g. about 10% for A, tt, K* , about 5% for the Kaons and 
nearly no change for the Ss. 
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FIG. 4: Effective masses for the different interacting chiral models and the ideal gas (vacuum 
values) case. The differences among the interacting approaches are less than 2%. 



These results, together with the steep contours from Fig. 2, suggest that the relative 
particle abundances "freeze" shortly after the spontaneous breaking of chiral symmetry. 
The success of our fit suggests extremely rapid chemical equilibration (through abundance- 
changing reactions) in the state with broken symmetry. Fig. 2 shows that the chemical 
composition of the hadronic system has to change substantially within a small temperature 
interval, just before freeze-out, even for the crossover transition (i.e. parameter set CI); for 
reference, we have indicated the dynamical path in the T — /i^ plane corresponding to the 
expansion of a perfect fluid (i.e. with constant entropy per net baryon |^4|)- While 2 — > n 
reactions are perhaps too slow to explain such rapid chemical equilibration [^, |3^ , m — > n 
processes with several particles in the initial state may be important as well ||3^, ^ |3^, 

Alternatively, the appearence of chemical equilibrium right after the phase transition (or 
the crossover) to the state of broken chiral symmetry might just be the outcome of the 
dynamical symmetry breaking process itself ETl, with statistical occupation of the various 



hadronic channels according to phase space ^ If so, number-changing reactions 

in the broken phase need not proceed at a high rate. To test this picture experimentally, 
it might be useful to consider central collisions of small ions like protons or deuterons, 
at similar energy and particle densities in the central region as for central Au+Au. For 
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systems of transverse extent comparable to the correlation lengths of the chiral condensates, 
the dynamical symmetry breaking process should be different from that in large systems 
(for example, the mean field approximation should not apply). The correlation lengths ^o-,c 
are given by 

, (Q/V) 

and accordingly for We evaluate the curvature of the thermodynamical potential at the 
global minimum and for T, fis, f^s at the freeze-out point. For parameter sets CI, CII, 
CIII wc obtain ^a- = 0.37 fm, 0.41 fm, 0.40 fm, respectively. For the correlation length of 
the strange condensate we obtain = 0.20 fm in all three cases. The correlation lengths 
are not very much smaller than, say, the radius of a proton. Thus, even if the freeze-out 
point for high-energy pp collisions happens to be close to that for Au-|-Au collisions at RHIC 
energies, the transition from the symmetry restored to the broken phase might be different. 
Finally, we also note that the correlation lengths obtained from our effective potential are 
not larger than the thermal correlation length 1/T at freeze-out, and so corrections beyond 
the mean-field approximation employed here should be analyzed in the future. 



CONCLUSION 



Particle ratios as calculated in a chiral SU {3)a — uj model are compared with RHIC data 
for Au-|-Au at ViSO AGeV and with noninteracting gas calculations. Since different versions 
of the chiral model show qualitatively different phase transition scenarios, we investigate 
whether the particle production, i.e. the chemistry of the system, is sensitive to the phase 
transition behaviour. Since we have shown that the current data are described by all three 
different phase transition scenarios and the ideal gas model, we can so far not favour or 
rule out any one scenario. In all interacting models the effective masses at freeze-out are 
shifted up to 15% from their vacuum values. The fitted chemical freeze-out temperatures and 
chemical potentials depend on the order of the phase transition. The crossover case yields 
15 MeV shifted T values as compared to the noninteracting gas model while the models with 
a first order phase transition yield more than 30 MeV lower temperatures. Furthermore, 
the fitted freeze-out points are located practically right on the phase transition boundary, in 
the first order phase transition scenarios, but T is always < Tg. This suggests that at RHIC 
the system emerges after the chiral chiral phase transition. This of course is only true if 
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a first order phase transition does actually occur in QCD at small chemical potentials and 
high T. "Freezing" of the relative abundances of various hadrons in the symmetric phase 
(at T > Tc) is excluded. 
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